Matrices, bases and matrix elements for cubic double crystallographic groups. 
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Matrices of the irreducible representations of double crystallographic point groups O, Td, O ® 
{1,/} and Td £8> {1,/} are derived. The characteristic polynomials (spinor bases) up to the sixth 
power are obtained. The method for the derivation of the general form of an arbitrary matrix 
element of a vector/tensor quantity is developed; as an application, the k ■ p matrix elements are 
calculated. It is demonstrated that the other known method for obtaining the bases of the irreducible 

representations of the double groups (L ■ S-diagonalization of a linear combination of spherical 
harmonics) is unreliable. 
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I. INTRODUCTION 

The studies of the part of the group theory which is 
used in condensed matter physics are commonly believed 
to be definitely finished by the middle of the XX cen- 
tury. Introduction to group theory is now an indispens- 
able part of condensed matter textbooks. Surprisingly, 
some important results (e.g. the matrices for the double 
irreducible representations) are missing. Consequently 
the power of group theory is not completely exploited. 

Condensed matter textbooks often contain character 
tables (without matrices) of the double groups, and this 
might give a wrong impression that a character table 
provides all symmetry group information one would ever 
need in condensed matter physics. This might be true for 
a non-degenerate case: using character tables one easily 
understands whether a matrix element (>I/|0|<I>) of some 
operator O between two non-degenerate states |W) and 
|<E>) must be zero by symmetry or not. The situation be- 
comes more complicated when the states \iff) and |$) are 
degenerate, and O is a vector (or a higher order tensor). 
In this case (^|0|$) is a matrix (or a set of matrices), 
and one would like to understand its structure rather 
than just revealing if it is zero or not. The knowledge of 
matrices allows one to get a deeper insight into (^|0|$): 
one can deduce its most general form using straightfor- 
ward formalism, see Sec. IIV Al 

Another usage of matrices is efficient constructing of 
the bases of irreducible representations (irreps). It is 
possible^ to construct a polynomial basis for an irre- 
ducible representation (irrep) without matrices: one finds 
a linear combination of the spherical harmonics (which 
correspond to the same degenerate level of a lonely atom) 

which diagonalizes the L ■ S operator (where L is the 
angular momentum, and S denotes spin). I see several 
drawbacks of this method though - unreliability, com- 
plexity, and narrowness: 

• Unreliability: I am not aware of a rigid proof of re- 
liability of this method (i.e., that every set of func- 
tions obtained using this method indeed does form 



a basis of an irreducible representation). To me 

the requirement of the L ■ S diagonalization seems 
too weak. My suspects are confirmed in Sec. IIV Bl 
where an example of a wrong basis obtained with 
this method is presented. 

• Complexity: In order to obtain a high (say, sixth) 
power polynomial basis of an irrep one has to di- 
agonalize a large matrix which may be hard to do 
analytically. 

• Narrowness: the method is restricted to polyno- 
mial functions. These can only approximate wave 
functions in the vicinity of zero, which is unsatis- 
factory, e.g., in a numerical calculation where one 
would like to reconstruct the wave function (or elec- 
tron density) in the entire atomic cell. The knowl- 
edge of matrices permits construction of projection 
operators^ which symmetry arbitrary (not nec- 
essary polynomial) wave functions. Such a sym- 
metrization would allow, for example, to get rid of 
unphysical parts of the wave function (which could 
appear, e.g., due to numerical errors), and might 
improve both accuracy and speed of numerical cal- 
culations. 

Recently Elder et ali derived general form of the k ■ p 
Hamiltonian for the case of T<j {1, 1} group (where 1 
and I are unity and spatial inversion operators, and <g> 
stands for the direct product) using physical approach 
(diagonalizing the Hamiltonian in the presence of spin- 
orbit interaction) . In this article I solve the same problem 
using projection operators approach^ which seems me 
simpler, easier to check by the reader, and more reliable. 

The calculation of matrices of irreducible representa- 
tions is based on the work of Dixoni, who is using the 
Burniside theorem, according to which every irreducible 
representation (irrep) of a group is contained in some di- 
rect product of certain number of its faithful (but not 
necessary irreducible) representations. Dixon^ demon- 
strated that in order to simultaneously block-diagonalize 
the set of direct products of matrices, it is enough to 
diagonalize one (specially prepared) matrix. 



2 



The practical realization of the Dixon's method, how- 
ever, may be problematic if the dimension of the direct 
product of representations is not small enough: the prob- 
lem arises due to the fact that it may be hard to diagonal- 
ize a large matrix analytically. I faced this problem when 
dealing with double groups with inversion, 0<g>{l, 1} and 
Td (g> { 1 , 1} , which have no faithful irreps in the standard 
(Pauli) gauge (which assumes that the inversion oper- 
ator I multiplies a spinor by —1). In order to resolve 
it I had to use somewhat less known Cartan gauge, see 
Sec. HIl This allowed me to obtain the matrices of the 
irreps analytically, but changed the double groups (in- 
cluding their character tables, cf. Table Q] and Ref. [f|). 

The choice of the gauge affects matrices of irreps and 
the character table but does not affect— physical quan- 
tities, e.g., bases of the irreps as well as matrix elements 
of operators (such as k-p) between the bands in a crystal. 
(Note, however, that matrices of an irrep are always de- 
fined up to a unitary transformation; so the correspond- 
ing basis sets are not unique.) 

In order to save space I have included only two gen- 
erators for each faithful irreducible representation in the 
text. However, the reader is encouraged to use supple- 
mentary material^, where matrices for all irreps sepa- 
rated in classes together with the transformation param- 
eters (see Sec. are available. 

It is important to check the obtained matrices. For- 
tunately, the check is much easier than the derivation; 
for convenience I provided a small program which cal- 
culates characters and multiplication tables of the irre- 
ducible representations. The program can be easily ex- 
panded by the readers for additional tests. It is available 
in the supplementary material^; there are versions for 
both maxima^ and Mathematica computer algebra soft- 
ware (CAS) systems. 



II. ISOMORPHISM BETWEEN CUBIC (O) AND 
TETRAHEDRAL (T d ) GROUPS 

The simple (geometrical) group O consists of all proper 
rotations^ of a cube and contains 24 elements divided in 
five classes as follows: 



class number — > 


1 


2 


3 


4 


5 


# of elements — > 


1 


3 


6 


6 


8 


rotation axis — > 




[110] 

a 


[100] 

a 


[100] 

a 


[111] 




m 


rotation angle — > 





7T 


tt/2 


7T 


2tt/3 


r$ character—)- 


3 


-1 


1 


-1 






(1) 



The last line in the Table ([TJ is the character of the (faith- 
ful) irrep, which is given by the usual 3D proper ro- 
tation matrices ([3]) which transform a cube into itself. 

The full tetrahedron symmetry group Td contains both 
proper and improper rotations [the latter are emphasized 
by the over line and red color in Table |(2J]. The corre- 



sponding 3D rotation matrices compose (faithful) T4 ir- 
rep. The group Td is classified similarly to the O-group: 



class number — > 
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2 


3 


4 
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II £ 1 j_ 

# of elements — s- 
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3 


6 


6 
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rotation axis — > 




[001] 

4 


[001] 

4 


[110] 

4 


[mi 

4 


rotation angle — > 





7T 


tt/2 


7T 


2tt/3 


Yi character— ^ 


3 


-1 


-1 


1 






where the last line is the character of the Y± irrep; dif- 
ferently from the Tr, irrep, T4 contains matrices having 
determinant —1. One can directly check that that the 
transformations of two simple groups O and Td obey the 
same multiplication table so that the two groups are iso- 
morphic. 

The matrices of the other irreps (ri,r2, and T3) are 
known^; otherwise one could deduce them from the direct 
products of the known faithful irreps using the Dixon 
method^. 

While the spatial inversion transformation operator / 
is uniquely defined in a 3D vector space (I = — 1), the 
transformation of a spinor might be different^— I mul- 
tiplies a spinor by a constant which may have any out 
of four values: ±1 (so called Pauli gauge) or ±i (Car- 
tan gauge). For the groups without inversion center [O 
and Td), the choice —1 is especially convenient since it 
is the only one which leads to isomorphism between the 
double groups^ Td and O. On the other hand, Car- 
tan gauge appears more convenient for studying double 
groups O ® {1,-0 and Td <8> {1,7} because it provides 
several faithful irreps to both groups. (In Pauli gauge 
these double groups posses no faithful irreps.) While 
I would prefer to perform all calculations in the (stan- 
dard) Pauli gauge, the advantage of Cartan gauge is cru- 
cial: it strongly reduces the dimension of the matrices 
which I had to diagonalize analytically while using the 
Dixon method^. For this reason, both gauges are used 
in this article: the results for double group O are ob- 
tained in Pauli gauge; the double group Td has been 
studied in both gauges: while matrices of the irreps and 
the character tables are gauge-dependent, the physical 
results arei^ gauge- invariant. Finally, I was unable^ to 
obtain analytical results for the double groups O £§>{!, 7} 
and Td <£> {1, 1} in Pauli gauge; for these last two double 
groups Cartan gauge (I = i in spinor space) was used. 

III. OBTAINING THE MATRICES FOR THE 
IRREPS 

An arbitrary element (transformation) of a double 
group can be characterized by the following (total five) 
parameters: (i) three Euler angles (a,/3, 7) which deter- 
mine the rotation ((3|) in the 3D space, (ii) the sign (or the 
branch number) in ([5]), and (iii) the presence/absence of 
inversion. 
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Below I provide a couple of generators for every faith- 
ful double irrcp. Complete sets of matrices grouped in 
classes can be found in the supplementary material^. 



A. Double groups O and Td 

In this section I use Pauli gauge so that the considered 
double groups are isomorphic and their matrices are the 
same (or similar). The matrices for the first five irreps 
are the same as for the corresponding single (geometric) 
groups, see, e.g., Ref. @- 

There are several definitions of Euler angles; I use the 
following connection between Euler angles and 3D rota- 
tion matrix, see ([111 1.4.54) or (01 -4.63): 



R 



3D 



R z 3B (a)R y 



(/3)^d(7), 



(3) 



where 




and 



(4) 



It is straightforward to obtain Euler angles for all proper 
rotations which transform a cube into itself. Substituted 
in ([3]), these Euler angles would produce 24 matrices of 
the IVrepresentation. Similarly, the set of matrices for 
the T4-irrep can be obtained from all rotations which 
transform a tetrahedron into itself. [If the rotation is 
improper, the corresponding 3D rotation matrix in ([3]) 
changes its sign.] 

Matrices for the Tq representation can be obtained by 
substituting the values of a, (3, 7 for the O-group into the 
expression ([l5|2.5.32) for the spinor rotation operator: 



B 1/J (a lJ 9,7) 



cos f e 4 
- sin I e l 



sin ^ e 1 
cos ^ e l " 



(5) 



where I have inserted =F which stands for two branches of 
the matrix function D 1 / 2 . [In physics textbooks a double 
group is often defined with the concept of "non-identical 
rotation by 27r" ; I prefer the more formal "two branches" 
definition-^ instead.] Note that in Pauli gauge used in 
this section the rotation parameters of the Td group gen- 
erate the same set of matrices, as the parameters for the 
O-group. 

I had to insert =p instead of ± in ([5]) in order to achieve 
the compatibility with Ref. @, where the class {6S4} in 
the character table of the double T^-group is understood 
to be composed of the ^irst-branch matrices. 

Both branches of (O produce 48 different 2x2 matrices 
which compose the so-called Tq representation. Using 
CAS one easily separates these matrices into eight classes 
and obtains the multiplication table which demonstrates 
that Tq is a faithful irrep. 



All irreps of a finite group are contained^ in a certain 
direct product of its faithful representations; this means 
that it is enough to know only one faithful (but not neces- 
sary irreducible) representation of a finite group in order 
to derive (at least numerically) matrices for all its irre- 
ducible representations. The search of irreps is further 
simplified if we know the character tabled for the double 
group, which tells us that all the missing (non-trivial) ir- 
reps (that is, T2, T3, Tr, Tg) are contained in the following 
direct products: 



r fi = r 7 



r 8 , r 4 ®r 5 = r 2 



r 4 



(6) 



The extraction algorithm (that is, the simultaneous 
transformation of all matrices of a reducible represen- 
tation into block diagonal structure) is invented and ex- 
plained by Dixoni. 

Since the double groups O and Td are isomorphic their 
irreps-matrices can be chosen to be the same. However, 
these two double groups are not identical (e.g., their ir- 
reps have different bases). The reason for this discrep- 
ancy is that coordinate functions are transformed differ- 
ently in these groups (for O, the corresponding irrep is 
while for Td it is T4). Due to the same reason the 
double groups O ® {1,1} and T d <£> {1, 1} also have dif- 
ferent bases, no matter what gauge (Pauli or Cartan) is 
used for the spatial inversion operator /. 

The generators of the faithful irreps Tq, T7, and T$ are: 



for r fi 



J_/l -1 

M 1 1 



and —= ( — 

V2 V 



(7) 



where \f±i = cxp[±i7r/4], 

• for T 7 : 

71 ( = 1 -l) 71 ( 

• for Ts: 



i V i 



(8) 




and 



(9) 



where every pair of matrices has the following transfor- 
mation parameters: 

1. a = 0,/3 = 7r/2,7 = 0, second (positive) branch 
in (O, and 

2. a = (3 = 7r/2,7 = it, first (negative) branch in <j5j> . 



4 



Double groups O <g> {1, /} and T d <g> {1, 1} 



for r 



All results of this section are obtained in Cartan gauge, 
see Sec.HU Both double groups, 0®{lj} and T d <E){lJ} 
have 16 classes (and the same number of irreps). The 
groups are not isomorphic - they have different char- 
acter tables (see Tab. [II] and Tab. IIII[) and non-similar 



matrices. The matrices for ten "single" irreps . . . 
can boH easily derived from the corresponding irreps 
. T 5 (the ones for the groups O and Td) as follows. 

. Tg", both first eight and 



For even representations^ 7 - 
last eight classes are given by the same matrices as for the 
group O, see Sec. IIII Al The same is valid for the "odd" 
representations r^...r^", except for that the matrices 
for the last eight classes are multiplied by —1. These 
"single" irreps arc identical for both double groups (in 
particular, the corresponding matrices are the same (or 
similar) for O ® {1, /} and T d <g> {1, /}). 

Without knowing apriori the character table generat- 
ing the irreps-matrices is somewhat more complicated. In 
case of O £g> { 1 , / } we depart from the transformation pa- 
rameters of the double group O; substituted in Eqs. ([3]) 
and (O, these parameters produce 48 (out of total 96) 
matrices of the irreps and Tg . In case of , the 
rest of the matrices is obtained by multiplying the first 
48 ones by —1; in case of r,l~, the multiplication constant 
is —i. 

Recursively going over different direct products of the 
known irreps we inevitably obtain the matrices^ for all 
16 irreps. The traces of the matrices composing the first 
8 classes coincide with one of the lines in Table H] the 
coinciding line determines the number (denoted by the 
subscript) of the newly obtained irrep. In addition to 
the number of the irrcp, we have to determine the par- 
ity (denoted by the ± superscript); I do this by build- 
ing polynomial bases and checking their parity under the 
transformation [x, y,z) — > (—x,—y,—z). 

The generators of the six faithful irreps (r^, , and 
Tg ) of the double group O (g> {1, 1} in Cartan gauge are: 



• for T+: 

J_/l -1 
V2 V i 

• for IV: 

Wl 1 

M- 1 1 

• for T+: 
- 1 



and 



and — = [ 

V2 V 




-\fi s/i 



and — = I 

y/2 V 



(10) 



(11) 



(12) 



2^2 




and 



(14) 



• for r 



8 



-2 



1 



1 



2^2 



V3\ 



2^/2 



2 i 
-1 iV3 1 \/3 
V\/3 i V3 -1 / 

/ \fi \J — 3i —\fi —y3i\ 
\fi V '— 3i —\f—i 
2-J^l -y/~i J^ii 



and 



(15) 



V 2VI -y/=zi -V^iJ 



where every pair of matrices has the following transfor- 
mation parameters: 

1. a = 0,/3 = 7r/2,7 = 0, inversion is absent, second 
(positive) branch in ([5]), and 

2. a = P = 7r/2,7 = 7r, inversion is present, first 
(negative) branch in ([5]). 

The corresponding characters are given in Table [TTJ 

The same procedure is used for calculating the irreps- 
matrices for the double group Td ® {1, /}, and the char- 
acters are given in Tabic IIIII The generators of the six 
faithful irreps (rjjS Yf, and Tf) of the double group 
T d ®{l,I} are: 



for r 



6 



-1 1 
-1 -1 



V2 

• for r 6 : 

i f-1 -1 

TlU -1 

• for T+: 

J_ ( 1 1 

V2 V- 1 1 



and * ^ 
y/2 V -V* 



and '(-£ ^ 

V2 V 



aild y/2 ( Vi Vi 



(16) 



(17) 



(18) 



• for r 

1 

71 



7 • 

-1 -1 



and 



1 

71 



i y/% 

1 \T—% 



(13) 



• for : 

i f 1 1 

71 V- 1 1 



and 



y/2 



yfi yfi 



(19) 



■5 



for 



2^2 




and 



(20) 



for r; 



l-2i 



2V2 






i — 


2i 


-\/3 -1 


V3 


2 i-\/3 



and 



\i\/3 -1 iV3 -i j 



( VI 



2V2 



— 3i —\i — v3i\ 

— 3i \A \/— 3i —\f—i 

2yf = l — 

V 2^ -V^3i — 



(21) 



where every pair of matrices has the following transfor- 
mation parameters: 

1. a = 0,/3 = 7r/2,7 = 0, inversion is absent, second 
(positive) branch in ([5]), and 

2. a = (3 = 7r/2,7 = 7r, inversion is present, first 
(negative) branch in ([5]). 



IV. GENERALIZED SELECTION RULES 
A. The calculation method 

Let us study the system which symmetry is given by 
some double group G which contains \G\ elements. It is 
convenient to present matrix elements ($|0|$) between 
degenerate bands in the matrix form. [For example, a 
matrix element between a (two-fold) I^-band and a (four- 
fold) I^-band is a 2x4 matrix.] Symmetry enforces re- 
strictions on matrix elements (\P|0|<I>) of an operator O 
between the bands and |$). Suppose that the states 
in the band ^ are transformed according to some irrep 
named A, the operator O is transformed according to 
some irrep B, and the states in the band are trans- 
formed according to some irrep C. 

The case when O is a scalar is trivial and will not be 
considered. In case when O is a vector— operator (e.g., 

O = O = p) the quantity (5 , |0|$) is characterized by 
three indexes 

(ijk) = (tfildjISfc) = J *t(A)6j# fc (A)dA, (22) 

where A represents all arguments of a wave function 
(except for spin) in some representation (e.g. A may 



be coordinate or momentum). In total (^E , |0| ( E > ) has 
n = dim A ■ dim B ■ dim C elements. The symmetry (i) re- 
quires that some of the matrix elements in (|22[) are zero, 
while among the others only few (sometimes - only one) 
are independent - see, e.g., Eq. (|34|) below. 

Under the action of a symmetry element g £ G the 
integrand in (|2"2"|) is transformed^ according to the direct 
product of three representations A ® B ® C: 

g(ijk) = £ A:„ l (g)B r3 (g)C k>k {g)(i'j'k'). (23) 

i',j',k' 



The group averaging operator 



sec 



(24) 



commutes with the integration in (|22[) and leaves all ma- 
trix elements intact: 



¥ G (ijk) = (ijk). 



(25) 



Let us associate every set of indexes (i,j,k) with some 
(ort honormal) bases element e/ of the n-dimensional com- 
plex space C": 



ej-B- (i,j,k) -> (ijk), l(i,j,k) = 
= (i - 1) • dimB • dimC + (j - 1) • dimC + k, 



(26) 



where we have taken into account the fact that different 
index sets (i, j, k) may correspond to the same values of 
matrix elements (ijk). 

The transformation (|23[) corresponds to a linear oper- 
ator in C n : 



g(ijk)= £ A* Vi (9)BrMCk>k(9)-(i'j'k') 

i',j',k' 



(27) 



where we used one-to-one correspondence between I 
and (i,j,k) defined in (|26|) . Similarly, the group aver- 
aging operator from ((24)) can be associated with a linear 
operator in C" : 

Pce^XyiPGlOe.', (/'|P G |Z) = j^-Y^ 1 ' 

l' 1 1 S 6G 

(28) 

Then we consider as a linear operator which acts in 
the space C": 



(29) 



A direct product of irreps is often reducible. Rcducibility 
means that there exists a unitary transformation which 
converts all (Vp G G) matrices (l'\g\l) from (|28p into the 
block diagonal form, and every block would correspond 
to some irrep of the group G. The group theory tells 



G 



us that the averaging | C | 1 X^eG^'I^IO destroys (av- 
erages to zero) all matrix blocks except for those (one- 
dimensional) ones, which correspond to the trivial irrep 
Ti (or r^). In other words, n-dimcnsional matrix Pq 
in (|29j) is similar to a diagonal matrix where the only 
non-zero elements are ones. 

This means that the operator : R™ — > R m is a 
projector, and that it has only two different eigenvalues: 
and 1. The degeneracy m of the eigenvalue 1 is the 
same as the number of times which the trivial irrep Ti 
(or r^~) enters in the direct product 



A <g> B Cg) C = mTi + other irrcps. 



(30) 



In case when m = 1, all matrix elements (ijk) are pro- 
portional to only one constant. The proportionality coef- 
ficients are just components of the eigenvector of the pro- 
jector Pg which corresponds to the eigenvalue 1. In or- 
der to obtain the eigenvector, we choose any set (i', j', k') 
which is not projected to zero by Pq. In other words 



(ijk) A* H {g)B j , :j {g)C k , k {g) 



(31) 



g&G 



where I can take any set of indexes (i',j', k') which pro- 
duces a non-zero value in pip . [In case when m = 0, any 
set (i',j',k') produces zero in (f3Tj) .] An example of the 
case m = 1 is the matrix element of k ■ p between the 
Tg conduction band and the valence band in GaAs. 
Since the symmetry of GaAs is described by the double 
group Td, the vector operator p is transformed according 
to the irrep T4. From the fact that 



r 7 = Tj. + r 3 + r 4 + 2r 5 



(32) 



we conclude that m — 1. One of the indexes set which 
produces a non-zero value in (|3"Tj) is k') = (1, 1, 1). 

Substituting this and two other index sets into (|31l) we 
obtain 

for (i',j',k') = (1,1,1) (ijk)<xi(k-ff)a 2 /6, 

for (i',j',k') = (1,1,2) (ijk)<x0, (33) 

for (i',j',k') = (1,2,1) (ijk) oc (k-ff)a 2 /6, 

where a = (<T2, c 2 , 03) is the set of Pauli matrices. Going 
over all possible values of (ijk) we conclude that in a 
zinc-blende structure the matrix element (Te\k ■ p ^7) is 
parametrized by one constant, as predicted in 



(r 6 |fc -p|r 7 ) oc (k ■ d)a 2 . 



(34) 



One may note that according to the bases written in 
Ref. [2(J, (Fe\k ■ p |r 7 ) oc k ■ a and ask why there is an 
extra 02 in (|3"3")l and ([54]) . This happens because my r> 
basis is different from the one in Winkler's book2£: 
in order to obtain Winkler's basis one has to (i) exchange 
basis functions in and (ii) change the sign in front of 
one of the basis functions. Since both of these operations 



correspond to unitary transformations of the irreps, both 
basis (pU)) and the matrix (j3"4"|) arc in agreement with the 
bases in Ref. poj . 

In case when m > 1, matrix elements (ijk) arc 
parametrized by m independent complex constants. I 
illustrate this case on the example of (Tg\k ■p\Ts). From 
the fact that 

r 8 g> r 4 ® r 8 = 2r a + 2r 2 + 4r 3 + er 4 + 6r 5 (35) 

we conclude that m = 2 in this case, so that 



(r 8 |fe -p|r 8 ) = c 1 m 1 + c 2 m 2 , cx,c 2 e 



(36) 



The generalized selection rules allow us to determine 
the matrices M\ and M 2 following the same prescription 
which we used above in order to obtain the result ([M)) 

for (Telk -p\Tj). Going over all possible values of (ijk) 
we obtain several (more than m = 2) linearly depen- 
dent^ non-zero matrices. Since we know that m = 2, it 
is enough to consider any two linearly independent ma- 
trices, e.g., For (i',j',k') = (1,1,2) we obtain 



Mi oc k\ 




(37) 



For (i',j', k') = (2, 1, 1) we obtain 



M 2 oc k\ 



V3\ 

-\/3 




(38) 



Any other set of values (i',j',k') produces^ a matrix 
which (together with Mi and M 2 ) forms a linearly de- 
pendent system of matrices. E.g., for (i',j',k') = (3,2,2) 



7 



class — > 


a 


6 


c 


d 


e 


/ 


5 


h 


-H- pit plpTnPTlt<5 V 


I 


5 


Q 


12 


8 


1 


5 


8 


rotation angle — > 







ir/2 


7T 


27r/3 





7r/2 


27r/3 


rotation axis — > 






A 




A 






A 


r. 


1 


1 


1 


1 


i 


1 


1 


1 




1 


1 


— 1 


— 1 


i 


1 


— 1 


1 


r 3 


2 


2 








-i 


2 





-1 


r 4 


3 


-1 


-1 


1 





3 


-1 





r 5 


3 


-1 


1 


-1 





3 


1 





r 6 


2 





-is/2 





1 


-2 


iV2 


-1 


r 7 


2 





iV2 





1 


-2 


-iV2 


-1 


r 8 


4 











-1 


-4 





1 



Table I. Character table for the double group Td in Cartan 
gauge. The classes are enumerated in the same order as in 
Ref. [f|. The characters for Ti . . . Ts and Fs are the same 
as in Ref. The characters for Te and ITV differ from the 
standard^ ones by an extra factor —i which appears in front 
of a/2- The 3D vectors are transformed according to the irrep 
r 4 ; for the spinors the appropriate irrep is Te. 



we obtain 



M 3 oc k\ 





f 


-V3i 





-3i \ 




-V3i 





3i 





1 





3i 









V -3i 





-V3 


i / 




( V 


\/3 


-3\ 









-3 


+k 2 





3 


V3 




V 3 


- 


73 0/ 


( 





2 V3? 













-2y/3i 


2v / 3i 











V 


- 


2\/3i 








(39) 



+k 3 



Any two matrix sets from {Mi, M 2 , M 3 } are linearly in- 
dependent, but the whole set {Mi,M 2 ,M 3 } is linearly 
dependent. This means, e.g., that the following expres- 
sion is equivalent to (|36j) : 



(T 8 \k-f\T 8 )=c' 1 M 1 



c' 2 M 3 , 



-It *-2 



c'n G 



(40) 



In case when the matrix element (Tg\k-p \T$) is calculated 
between the same I^-bands, the two complex constants 
in (j3"6")l or in (|4"U)) gain a constraint: their values should 
be chosen in such a way that (r§|fc ■ p \T%) is Hermitian. 



In sections IWTTT1 IVI B 21 fVI B 3) and lVIB4I I present 
matrix elements of the k ■ p operator which are ob- 
tained using CAS according to the method described in 
Sec. lTVAl 



B. Comparison with Elder et al*i 

Elder et alX recently found k ■ p matrix elements using 
Lowdin approach which allows to split (approximately) 
the Hamiltonian into several blocks in such a way, that 
each block corresponds to a separate energy level (or to 
an irrep). 

Let us compare the result (|36|) with ( 31 ) . We have 
to take into account the difference between the bases 
on p. [l| 19 and (JTSD. One notices that the linear basis 
from (|70[) is equivalent but not equal to the linear basis 
for the irrep T% in Ref. [l[ . From the bases comparison I 
conclude that in order to translate my Tg irrep into the 
notations of Ref. [l| I have to apply the following unitary 
transformation to my T$ matrices: 



r' = ur*u~ 



u = 








-1 
















1 











\o 


-1 









(41) 



With the matrices T' s I do recover the linear basis for T s 
on p. 019: 

[V3(iy + x) t,-2zt+(iy + a;) I, 

-2z I +(iy - x) t,V3(iy - x) J.]. 



(42) 



However, instead of ([l|30a) I obtain (in the notations of 
Ref. Q): 



K 



r 8 ,r 8 oc 



+const 



( 


3k+ 





-k- 





-A/3fc 4 


-4fc 3 


VZk + 





\V3k+ 





3fc„ 


( 


3fc+ 


6fc 3 


-5fc_ 





V3k+ 


-2*3 


\/3fc_ 





\-V3k+ 


-6fc 3 


3k- 



V3k-\ 
4fc 3 
-k+ 
/ 

-V3k-\ 
2k 3 
-5fc+ 
/ 



(43) 



which does not fully agree with (030a). 

Next, with the matrices T' 8 I obtain [apart from (|42[)] 
two second-order polynomial bases which are incompat- 
ible with the ones on p. [I|l9: 



and 



[(2z 2 -y 2 -x 2 )i,V3(-y 2 +x 2 )t, 
V3(y 2 - x 2 ) U-2z 2 + y 2 + x 2 ) t], 

[-V3(y + ix)z f, -{y + ix)z 1 +2xy f, 

(y — «ar)z f +2xy I, Vi(y — ix)z I}. 



(44) 



(45) 



The reason for this incompatibility might be the fact that 
the sixth basis on p.[l|l9 



ipi = -2yz t +3iy 2 I +2ixz f -Axy 1 -3ia; 2 1, 
^2 = i%/3(2z 2 t ~2%z i -y 2 t -2xz | -a; 2 f), 
(^ 3 = i\/3(-2z 2 I +2iyz t +y 2 I -2xz f +x 2 I), 
<y5 4 = 2yz 4I -3iy 2 t +2ixz 1 -4xy t +3i.t 2 f 



(46) 
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If 


i 


1 


1 


1 


1 


1 


1 


1 


±1 


±1 


±1 


±1 


±1 


±1 


±1 


±1 


It 

2 


i 


1 


-1 


-1 


1 


1 


-1 


1 


±1 


±1 


Tl 


Tl 


±1 


±1 


Tl 


±1 


if 


2 


2 








-1 


2 





-1 


±2 


±2 








Tl 


±2 





Tl 


if 


3 


-1 


-1 


1 





3 


-1 





±3 


Tl 


Tl 


±1 





±3 


Tl 





it 


3 


-1 


1 


-1 





3 


1 





±3 


Tl 


±1 


Tl 





±3 


±1 





it 


2 





y/2 





1 


-2 


-y/2 


-1 


±2i 










±i 


T2i 


W2i 


T* 


if 


2 





-V2 





1 


-2 


V2 


-1 


±2i 





Tv^z 





±i 


T2i 


±y/2i 


T« 


it 


4 











-1 


-4 





1 


±4i 











T'i 


T4i 





±i 



Table II. Characters of the irreps of the double O ® {1,1} group in Cartan gauge. (The inversion operator I multiplies a 
spinor by i.) The superscript ± in the irreps-notation means that all its polynomial bases are even/odd with respect to 
the transformation (x,y,z) — > (—x,—y,—z). The 3D vectors are transformed according to the irrep for the spinors the 
appropriate irrep is TJ. 



if 


1 


1 


1 


1 


1 


1 


1 


1 


±1 


±1 


±1 


±1 


±1 


±1 


±1 


±1 


if 


1 


1 


-1 


-1 


1 


1 


-1 


1 


±1 


±1 


Tl 


Tl 


±1 


±1 


Tl 


±1 


it 


2 


2 








-1 


2 





-1 


±2 


±2 








Tl 


±2 





Tl 


if 


3 


-1 


-1 


1 





3 


-1 





±3 


Tl 


Tl 


±1 





±3 


Tl 





it 


3 


-1 


1 


-1 





3 


1 





3 


Tl 


±1 


Tl 





±3 


±1 





1 6 


2 





-y/2i 





1 


-2 


V2i 


-1 


±2i 





T^2 





±i 


T2i 


±^2 


Ti 


1 7 


2 





y/2i 





1 


_2 


-V2z 


-1 


±2i 





±y/2 





±i 


T2i 


W2 


T« 


1 8 


4 











-1 


-4 





1 


±4i 











Ti 


T4i 





±i 



Table III. Characters of the irreps of the double Td <8 {1, 1} group in Cartan gauge. The 3D vectors are transformed according 
to the irrep rj; for the spinors the appropriate irrep is Tg". 



is wrong. (All other bases on p.[l|l9 are correct^) First 
let me express the basis (|4T>|) as a linear combination of 
spherical harmonics^: 



'— [2y 2il t + (y 2 , 2 



Y'2 



V?3 



5^2,-2) I] 

-4iy| (V2Y 2A I +\/3r 2 , t) 

5^2,2) t] 



(47) 




Next, one can checb^ that the system of functions (|47j) 

is (i) orthogonal^ and (ii) diagonalizes^ the L ■ S op- 
erator. According to the commonly accepted belief, this 
should mean that (jT7|) is a basis of an irrep of the two 



groups: O® {1,1} and T^® {1,1}. However, the direct 
check-i 3 - (see the file badBasis.txt in the supplementary 
material^) shows that the system of functions (jT7|) docs 
not form a basis of any of the two considered groups. In 
particular, the rotation around the j/-axis by tt/2 (which 
is a symmetry transformation of both groups) does not 
transform the system of functions (|47l) into linear combi- 
nations of themselves. 

I am not aware of a rigid proof of reliability of the stan- 
dard LS^diagonalization method for finding the basis of 
an irrep. Moreover, the system of functions (|47|) is an ex- 
ample which demonstrates unreliability of the commonly 
accepted method of obtaining the polynomial bases of 
irreps of double groups. 

On the contrary, the situation with the method of pro- 
jection operators^ (used in this article) is clear: pro- 
jection operators can never produce a wrong basis of an 
irrep. 



V. CONCLUSION 

From the application point of view, the most impor- 
tant result of the article is the general form of the matrix 
elements of the k ■ p-operator. I expect it to be useful for 
studying effects in many important materials, for exam- 
ple, in 



I 

• zinc-blende semiconductors (GaAs, InAs, etc.), 

• semiconductors with inversion center (Ge, Si, etc.), 
and 

• cubic perovskites (SrTi03, LaA103, etc.). 

There are, however, other important (and more funda- 
mental) results discussed in the next two paragraphs. 
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a. The matrices of the irreducible representations of 
the double cubic groups O, Td, 0®{1,I} and Td® {1, 1} 
are obtained. First, these matrices are important for 
the derivation of the bases of irreducible representations, 
which is useful both in analytical and in numerical cal- 
culations. This is the only reliable method for the bases 
derivation I am aware of. An alternative approach is used 
in Ref. [l| and is based on the statement (which seems to 
be commonly believed) that any linear combination of 
spherical harmonics Yi m with the same value of I which 

diagonalizcs the i-5-opcrator should be a basis for some 
irreducible representation. In Sec. IIVBI I disprove this 
statement using a counter-example. There is somewhat 
similar method used in Ref. [U to construct a linear basis 
for the Tg-irrep of the double group Td, which assumes, 
however, that at least one basis function is known or can 
be guessed. Having tried to guess elements of the poly- 
nomial bases which I have obtained (see Sec. IVI CI and 
Ref. Q) I conclude that this "guessing" requirement is 
almost never fulfilled, so the method^ is practically use- 
less. 

b. An algorithm for obtaining general matrix struc- 
ture of matrix elements between degenerate energy lev- 
els is developed. This algorithm is more reliable than 
obtaining the matrix elements from the some manually 
(casually) chosen set of basis functions: there is always 
a chance, that due to the oversimplified (not general 
enough) basis one obtains wrong matrix elements: e.g., 
considering linear basis for the Ts-irrep of Td group, one 
might erroneously conclude that (rg|fc • p \ Y%) =0. 

This algorithm might be especially useful for analytical 
calculations with the k ■ p method; for this purpose the 
general form of the k ■ p matrix elements is derived (see 
Sec lVIBl in the Appendix). 

Note that there is an inconsistency between the results 
of this article and Ref. [2(| : the 8x8 Kane model in Ta- 
ble [2(|C.8 which has too many parameters (while only 
two are allowed) . There is also an inconsistency (of both 
my results and Table [2(|C.8) with matrix elements ob- 

I 



tained in Ref. but I believe that it is connected with 
the invalidity of the basisi, see Sec. IIVBI 

As a final note, let me mention somewhat misleading 
notation Oh which is often used for both O ® {1, /} and 
Td® {1,1} double groups. On the first glance it makes 
no sense to consider the double groups O ® {1,1} and 
Td <8> {1, /} separately: for example, in Pauli gauge these 
groups are isomorphic; moreover, any polynomial basis of 
some irrep of the double group O ® {1,1} is also a basis 
of some irrep of the double group Td® {1,1}. However, 
these bases often correspond to different irreducible rep- 
resentations, so it would be incorrect to claim that two 
double groups O ® {1,1} and Td ® are identical 

(no matter what gauge is used). See also the note in the 
end of Sec. IVI C 21 

VI. APPENDIX 
A. Notations 

I assume that \f±i = exp[±i7r/4], (01, a-i, a 3) are the 
usual Pauli matrices, N is the set of all positive integers, 
and C is the set of complex numbers. By saying that "two 
matrices are similar" I mean that they are connected by 
some similarity transformation. I denote 1 as an identity 
operator, and / as the spatial inversion operator. 

The irreducible representations are named according 
to the following rules: 

• An irrep named Tq (or Tq for the groups with in- 
version) should have a basis [f,i]- 

• All polynomial bases of "even" irreps (marked with 
the "plus" sign) must be invariant with respect to 
the transformation (x, y, z) — ¥ (—x, —y, —z). 

Equations from external sources are cited, as 
([citation]NN), where NN is the equation num- 
ber. For example, ( [l| 31) stands for the "equation (31) 
in the article fif" . 



B. General form of the k ■ p matrix elements 



1. k ■ p matrix elements for the double group O 



<r 6 |fc-p|r 6 ) oc k-a, (T 6 \k-p\r 7 ) = 0, (r 7 \k ■ p\r 7 ) oc k ■ a* , 



, V3 1 0\ f-V3i i \ (0 2 



(48) 
(49) 



V3 1 



-n/3 



-2 (T 
002/' 



(50) 
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(r 8 |A; -p\T s ) cx fci 



/ V3 -1\ / -y/3i 

V3 1 

1 \/3 

\-l \/3 / 



+ fc 2 



-const 



fci 



-\/3 0-3 
0-3 -V3 
V -1 -V3 j 



V5i 





1 \ 




n o 








U 


—i 


U 




n o 
U z 


U 


u 


j 


o 







-2 





u 








vo 


o 


-2 


V3i 





-i \ 









°\ 





3 / 





+ h 









-3i 







° 





=) 





-V3i 


) 




V o 








+ 



(51) 



2. fc • p matrix elements for the double group Td 

I have obtained bases and matrix elements of the double group Td in both Pauli and Cartan gauges (see Table U for 
the Cartan characters). The results are gauge- independent, as predicted by theory 



<r 6 |fc-p|r 6 > = o, (T 7 \k-p\r 7 ) = o. 

Matrix elements (MEs) (r 6 |fc ■ p\T 7 ) are written in (|34|) . 



(r 6 |fc -p\T s ) cx ki 



1 x/3 
-V3 -1 



-i V3i 
-i 





















o 



(52) 



(53) 



(tVlMr.>ocfa(_°5j; f 

Matrix elements (r§|fc • p |r§) are written in (|36l) . 



i V3i 
V3 i -i 



0-20 
2 



(54) 



3. k ■ p matrix elements for the double group O eg) {1, /} 
Matrix elements between both odd or both even parity states are zero, e.g., 

(r+|fc.^r+) = (r 6 -|fc.pir7) = o. 

The other matrix elements are 

(T%\k -p\Tj ) = 0, (r£\k-p\Ve) cx (k ■ a)ai, (rf \k ■ p\Tj ) cx (kiai + k 2 a 2 - ^30-3)0-2, 

, r+| r . A/3 -i 0\ , , /0 2 \ / V3 -1 



(55) 

(56) 
(57) 



/r+if -ir-\ j. f 1 iV3 0V , ^ 
(rf |fc-p|r 8 ) a fei [ „ „ „„I+K2 



2 



-i iVi) ^ fc3 \\ -iy/3 



1 V3 



(58) 



const • 



fci 



/V3 


/ 





\ 




f° 


-2 





°) 




' V 








- 2 \ 








V3 


-1 







- 


-t-\/3 


—i 




fc 3 


-V3 


i 














/ 


iV3 


+ fc 2 








-1 


V3 


■+ 


— i 


V3 








V-i 


-V3 





/ 




^ 








oy 




^ 





2i 




-2 





<M 




/0 2i 









/ 


- 


v^i 










V3i 
-1 


i 

V3 


+ fc 2 










-V3 
—i 


-1 

V3i 


+ &3 




-2 


-2 













2 









\2i 








J 









1 







(59) 
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4- k ■ p matrix elements for the double group Ta ® {1, /} 
Like in Sec. IVI B 3[ many MEs are zero by parity, and we do not write them here. The other MEs are 

(r+|£.|-|i7) = o, (r+|fe.#|rf> = o, 

(r+|fc • p|rf) cx (k ■ a)crx (r^|fc -p\T^) cx (fcidi + k 2 a 2 - k 3 a 3 )a 2 , 



(60) 
(61) 



/r+iiT -lr-\ fc ^ -i 0\ , , /0 2 \ / V3 -1 



(62) 



1 iV3 
2 



{T+\k .j?|r 8 ->ocM;; 'n" o n)+M ~ q ^ ^J+^U _^ g q 



2i 



1 v/3 



(63) 



/\/3 



(r+ifc.pir-)^^ 



o o \ 



-const 



\/3 -1 
i 
\-i -\/3 / 

^ -2 ^ 

V3i i 
-1 V3 

\-2 oo o y 



+ fc 2 



A) -2 \ 

-iV3 -i 

-1 V3 

\2 o o o y 



/ -2\ 



-V3 i 

i V3 

\ o o 2z o y 

/0 2i \ / -VZi -» \ 
-V3 — 1 , 0-2 
-i V3i 3 -2 

\2i / \ 1 -\/3/ 



(64) 



C. The polynomial bases 



With the information from Sec. Ill Al it is straightforward to get the bases for the irrcps using the standard projection 
operators technique (see, e.g, pp. ^82-83 or pp.Q64-65). 

First, let me mention some (scalar) characteristic polynomials of the irrep T±: 

^1 = 1, Vi 1 = x 2n + y 2n + z 2n , ^ 2 =y 2n z 2 + y 2 z 2n + x 2n z 2 + x 2n y 2 + x 2 z 2n + x 2 y 2n , ^\ = (V> 2 ) 2 ", (Tj) 

where n £ N is a positive integer—, and t/j 2 denotes an arbitrary characteristic polynomial of the I^-irrep, e.g., 

/ 2 ; 2 2 2n 2 2n 2n 2 , 2n 2 , 2 2n 2n 2 /"pi \ 

V'J = a^yz, ^2 = V z — x z — y z + x z + x y — x y . (1 2 ) 



Every basis below remains valid if all its components are multiplied by any expression from jr/iP . The presence of 
such bases is considered as obvious, so they are not explicitly mentioned below. 

In the main text of the article I write only few simplest polynomial bases for each irrep. (In the supplementary 
material^ the polynomial bases are available up to the sixth order.) 

Three-dimensional vectors are transformed differently in O and Td groups. In the O-case, the corresponding irrep 
is T^; in the T^-case, 3D-vectors are transformed according to the irrep IV (See the captions of the Tables HT1 and Hill 
for the groups with inversion centers.) As a result, despite the fact that the two groups are isomorphic in Pauli gauge, 
they still have different bases and different k ■ p matrix elements. 



1. Polynomial bases of the double group O 



For T 6 : 



[1,1], [-z 2n - x t -iy 2 ^ 1 I -x 2 "- 1 l,^ 2 ' 1 - 1 I +iy 2n - 1 t -x 2n ~ x % n e N, 

4} 2 {yz l +ixz I +xy f, yz j" — ixz j" —xy 4-], 



(65) 
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where ip 2 denotes an arbitrary characteristic polynomial of the I^-irrep, taken, e.g., from fll^l ). For Tj: 

[-(yz) 2 "" 1 t +i(xz) 2n - 1 t + (xy) 2 "- 1 |, (yz) 2 ™- 1 | +z(xz) 2 ™- 1 | + (xy) 2 "- 1 fl. 



(66) 



For Ts 



[-^(x 2 ™- 1 + ly 211 - 1 ) t^z 2 "- 1 1 -iy 2 ™" 1 1 -x 2 " 1 - 1 |, 

2Z 2 "" 1 | -iy 2 ™" 1 t +x 2n - 1 t, V^z 2 " -1 - iy 2 ™" 1 ) ]], 
ip 2 [2z I +x t -iy t,V3(-x + iy)l, V3(x + iy) f, 2z t -a; | -iy |] , (67) 

ip 2 [V?>z{y + ix) *f,yz I +ixz I — 2xy f, — yz j" +ixz j" — 2xy \/3z(— y + ix) \], 
[V3(-y 2 ™ + x 2n ) \.,(-2z 2n + y 2n + x 2n ) f, (2z 2 ™ - y 2n - x 2n ) |, V3(y 2n - x 2n ) f]. 

U. Polynomial bases of the double group Td 



For r 6 : 



[t, {] , V- 2 [t, I] , V> 2 [z 3 t +V 3 t +* 3 t, z 3 I +y 3 I +x 3 |] , 

[yz 2 ™ | + iy 2n z t +ixz 2n | -ixy 2 ™ | -ix 2 "z f -x 2 "y |, (68) 
- yz 2n t -iy 2 "z | +ixz 2 ™ t -ixy 2 ™ t +ix 2n z | +x 2 "y t]- 
For r 7 : 

[z 2 ™- 1 | + i y 2 "- 1 t -x 2 "- 1 t, z 2 "- 1 t +*y 2 "" 1 I +-T 2 "- 1 |], [-yz t +ixz t +xy |, yz | +*xz | +xy t], 
^[z 2 "- 1 | +J), 2 ' 1 - 1 | -a; 2 "- 1 t, z 2 "" 1 t +Jy 2 "' I +i 2 " -1 I], Tp 2 [-yzt+ixzt+xyl,yzl+ixzl+xy f]. 
For T 



[-yz t +ixz f -2xy |, \/3z(y - x) I, - V3z(y - ix) f, yz | +ixz | -2xy t], 

V-> 2 [— yz j" +ixz j" — 2xy J,, V3z(y — x) — Viz(y — ix) f , yz +ixz J, —2xy f], 

[2Z 2 ™- 1 1 -iy 2 ™" 1 1 +X 2 ™- 1 t,\/3(iy 2n_1 - a; 2 ™" 1 ) |, 

^(izj 2 "- 1 +X 2 "- 1 ) t^z 2 ' 1 - 1 t -iy 2 ' 1-1 I -x 2 "- 1 |], (70) 

[V3(-y 2 ™ + x 2 ") |, (-2z 2 ™ + y 2n + x 2n ) f, (2z 2 ™ - y 2 " - x 2 ™) |, v^y 2 ™ - x 2 ™) f], 

[2yz 2 ™ t +iy 2n z | +2 Z xz 2 ™ | -ixy 2 " f +zz 2 ™z | - x 2 "y f, V3(-zy 2n z t -ixy 2 " | +ix 2 "z t -z 2 "y |), 

\/3(-iy 2 "z | +ixy 2 " t +ix 2n z I -x 2 "y t),2yz 2n | +iy 2n z t -2ixz 2n | +ixy 2n | +ia; 2n z t -a; 2n y |]. 

Note that the two bases (|6"Tj) and ([70)) are not equivalent, because the order of the functions in an irrep-basis matters: 
e.g., one could interchange elements in the first basis set in (|70| in such a way that it becomes equal to the first basis 
set in dHTJ; however, such an exchange would make second basis sets in ([70]) and in (jBT)) different. Generally, a wave 
function is given by a mixture of several basis sets. 

3. Polynomial bases of the double group O ® 

In case of O ® {1, 1} the vectors/spinors are transformed according to the /Tg irreps. The bases for Tq are: 
[t, |], [yz 2 ™- 1 | -y 2n - x z I -ixz 2n - x | +xy 2n - 1 t +ix 2n - 1 z | -x 2 ™"^ f, 
yz 2 ™' 1 t -y 2 ""^ t +ixz 2n - 1 t -xy 2 ™ -1 I -ix^^z f +x 2 "- x y |]. 



(71) 



For Tg 



[z 2 ™" 1 | +iy 2 ™~ 1 t -x 2 ™- 1 f, - z 2n ~ x t -iy 2 "" 1 I -x 2 ™- 1 |], 
[-yz 2 ™ t -ixz 2 " t -ixy 2 ™ f + ^ 2 ™^ I +ix 2 "z | -x 2 "y f, 



yz 2 " | -iy 2n z t -ixz 2 " | - Z xy 2 " | -ix 2 "z t +x 2 "y |], 
V> 2 [yz t — ixz t — xy |,yz +ixz J, +xy f]- 



(72) 



For YJ: 



For Tj 



For T+: 



[(yz) 2 "" 1 I +z(xz) 2 "- 1 | +(xy) 2 "- 1 f, - M 2 "" 1 t +^(xz) 2n - 1 t +(.Ty) 2 "- 1 |], 
[yz 2 "- 1 i +y 2n - 1 z i +ixz 2n ~ 1 | + xy 2 ^ 1 f +ia; 2Tl - 1 z | W^y t, 

-yz 2 ™- 1 t -y 2 ™^ t + ixz 2 "- 1 t | + Z a; 2 "- 1 z f +x 2 "^ 1 y |], 



i' 2 [I, t] . [ - 2/* 2n t -«y 2 ™2 | +ixz 2n t -i.Ty 2n t +ix 2 "z I +x 2 "y f. 

yz 2 " | +iy 2n z t +ixz 2 " | -ixy 2 " I ~ix 2n z f ~~x 2n y |]. 



[-z 2 " t +2y 2 " t -x 2 " t,z 2n I -2y 2n | +x 2n i,iV3(-z 2n + x 2n ) |,i%/3(z 2 " 1 -x 2 ") f], 
^. 2 [V3(z 2 "- 1 t ~x 2n ' x i)^{z 2n - x | +VZx 2n - 1 f), 

iz 2 "- 1 | +2y 2 "~ 1 t -iT 2 "- 1 t,^ 2 "- 1 t +2y 2 "~ 1 | +iT 2 '"- 1 |], 
[V3(yz 4- ~xy \), — Vs(yz j" +xy J,), «yz j" — 2.tz j" — ixy J,, — iyz 4- — 2.tz 4- — ixy f] 



For T 
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[VZ{z 2n - 1 I +X 2 "- 1 t),^ 2 "" 1 | +2y 2 "~ 1 t -ix 2 "- 1 t, 

^(z 2 "- 1 1 -iy 2 " -1 1), 1 -*y 2 ' 1-1 1 +2X 2 ' 1 - 1 I], 

V> 2 [z 2 | -2y 2 I +x 2 \.,iV3(-z 2 | +x 2 ) I, z 2 t +y 2 t -2z 2 t, V^fz 2 - y 2 ) t], 

ip 2 [V3(yz j" +a:y 4-), — iyz j" +2xz j" +ixy 4-, \/3(— iccz 4- +xy t), 2yz 4- — ixz 4- — xy f\, 
[V3(-yz 2n t -ixz 2n t - zx 2 "z 4, +x 2n y f), -^yz 2 " t +2y 2n z | +xz 2n f -2a;y 2 " t -z 2 "z I -ix 2 "y f, 
V3(-yz 2n 4, -iy 2n z f + ixz 2 " | -izy 2 " |), yz 2 " | -m/"z t -ixz 2n | -ixy 2 " | +2ia; 2 "z | -2.T 2 "y 1] 

^. Polynomial bases of the double group Td <8> {1, 7} 



In case of ® {1,7} the vectors/spinors are transformed according to the T 4 /Y^ irreps. The bases for Y 

[t, 44, [yz 2 ™' 1 I -y 2 "-^ I -ixz 2n - x I +xy 2 ™- 1 t +^ 2 ™" 1 z 4. -x 2n ^y f, 
yz 2 "- 1 t -y 2 " _1 z t +WZ 2 "- 1 t -a;y 2n_1 I -^"^z t +^ 2 ™ _1 y I]- 



For 



For T^: 



For Ff 



^ 2 [;,t], [- y^ 2 ™ t -*y 2 "z I +ixz 2n t -i.xy 2n t +ix 2n z I +x 2n y f, 
yz 2,i I +^y 2n z t +ixz 2 " | -ixy 2n | -ia; 2 "z t -s 2 "y ]]. 



[yz 4- +ixz 4- +a;y t, — y z t +«a;z j" +xy 4,], 

V- 2 ^ 2 ™" 1 1 +«y 2 ™ -1 1 +X 2 ™- 1 1, z 2 ™- 1 4. +iy 2n - 1 1 — .x 2n ^ f]- 



z 2 "- 1 4- 


+iy 2 ™" 1 t - 


-x 2 "- 1 1, - 


2 2n-l 


t 


-*y 2 - 


u 




[-y^ 2n 


t +*y 2 "^ 4- 


-ixz 2n t - 


ixy 2n 


t 


+ix 2n z 




-^ 2 ™y t, 






yz 2 " | - 


iy 2n z 


t 


2n 

-zxz 


4-- 


^xy 2n 4--zx 2n zt+^ 2 "y|], 



ip 2 [yz f — ixz f — sy JoJ/z 4- +ixz I +xy f]- 
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For T+: 



[-z 2n f +2y 2n | -x 2n t, z 2n 1 - 2y 2n | +x 2n |, iVi{-z 2n + x 2n ) |, iV3(z 2n - a; 2 ") f], 
^ 2 [v / 3(^ 2 ™- 1 t "X 2 ™- 1 l),^ 2 "- 1 | +X 2 "- 1 t), 

iz 2 ™" 1 | +2y 2 ™- 1 t -ix 2n - x t, ^ 2 " _1 t +2y 2n ~ 1 I +ix 2n ~ 1 I], (81) 
[^{(yz) 2 "- 1 I -(xy) 2 "- 1 t}, - ^{(yz) 2 "- 1 t Hxy) 2 ^ 1 |}, 

z(yz) 2 "- 1 t - 2(xz) 2 "- 1 t -tixy) 2 "- 1 |, ^{yz) 2 ^ 1 1 -2(xz) 2 "" 1 | -z(xy) 2 ' 1 - 1 t]. 



For To 



[Vi(z 2n - 1 I +x 2n ~ X t),jz 2 " -1 I +2y 2 "^ 1 t -ix 2 "- 1 t, 

^(z 2 "- 1 t -V^y 2 "" 1 |), -z 2 "- 1 t -«y 2 " -1 I +2Z 2 "- 1 J,], 
^ 2 [z 2 | -2y 2 | +x 2 |^V3(a; 2 - z 2 ) |, z 2 t +y 2 t -2x 2 f, \/3(z 2 - y 2 ) f], 

V> 2 [\/3(yz j" +xy 4), — iyz j" +2xz j" +ia;y 4, v3(— 4- +xy t), 2yz 4- — ixz 4 — xy f]. 



(82) 
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